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ABSTRACT 

We present wide-band radio observations spanning from 1.4 GHz to 350 GHz 
of the afterglow of GRB 991216, taken from 1 to 80 days after the burst. The 
optical and X-ray afterglow of this burst were fairly typical and are explained 
by a jet fireball. In contrast, the radio afterglow is unusual in two respects: 

(a) the radio light curve does not show the usual rise to maximum flux on 
timescales of weeks and instead appears to be declining already on day 1 and 

(b) the power law indices show significant steepening from the radio through the 
X-ray bands. We show that the standard fireball model, in which the afterglow 
is from a forward shock, is unable to account for (b) and we conclude that the 
bulk of the radio emission must arise from a different source. We consider two 
models, neither of which can be ruled out with the existing data. In the first 
(conventional) model, the early radio emission is attributed to emission from 
the reverse shock as in the case of GRB 990123. We predict that the prompt 
optical emission would have been as bright (or brighter) than 8th magnitude. In 
the second (exotic) model, the radio emission originates from the forward shock 
of an isotropically energetic fireball (10 54 erg) expanding into a tenuous medium 
(10~ 4 cm -3 ). The resulting fireball would remain relativistic for months and is 
potentially resolvable with VLBI techniques. Finally, we note that the near-IR 
bump of the afterglow is similar to that seen in GRB 971214 and no fireball 
model can explain this bump. 
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1. Introduction 

The intense gamma-ray burst GRB 991216 was detected on 1999 December 16.67 UT 
by the Burst and Transient Experiment (BATSE) on board the Compton Gamma Ray 
Observatory satellite ( Kippen, Prcccc fc Giblin 1999|) . Follow-up observations with the PCA 



instrument on the Rossi X-ray Timing Explorer (RXTE) satellite resulted in the detection 
of a previously uncataloged X-ray source, which was subsequently seen to fade by a factor 
of five, seven hours later ( Takcshima et al. 1999 ). Uglesich et al. (1999) identified a fading 



optical source, at a position consistent with the RXTE transient, and shortly thereafter the 
radio counterpart was discovered fllaylor fc Berger 1999| ). 



Here we present radio measurements of this burst from 1 GHz to 350 GHz. While 
the emission from X-ray and optical afterglow was fairly typical flHalpern et al. 2000|) , the 
radio afterglow of GRB 991216 was unusual in two respects. First, the onset of the decay 
began much earlier than that in most radio afterglows. Second, the temporal decay indices 
in the radio, optical and X-ray bands are markedly different from each other. We explore a 
number of possible explanations for these behaviors. 



2. Observations 

Very Large Array (VLAf\): A log of the observations and flux density measurements are 
summarized in Table 0. We used J0509+1011 (at 8.46 GHz and 4.86 GHz) and J0530+135 
(at 1.43 GHz) for phase calibration. J0542+498 was used for flux calibration at all 
frequencies. 

Very Long Baseline Array (VLBA 6 ): A single 2-hr observation was carried out at 8.42 GHz 
and 2-bit samples of a 64 MHz bandwidth signal in one hand of polarization were recorded. 
The nearby (< 1.1°) calibrator J0509+1011, a core-jet source, was observed every 3 minutes 
for delay, fringe rate and fringe phase calibration. The total flux density of the calibrator 
was found to be 9.5% less than was measured by the VLA on the same day. Given that the 
jet of J0509+1011 is likely to have some extended emission that is not detectable by the 
VLBA, it is likely that the absolute flux calibration of the VLBA is well within its nominal 
value of 5%. 

The radio afterglow was detected at a position of (epoch J2000) a = 5 h 9 m 31. 2983 s , 



6 The NRAO is a facility of the National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. NRAO operates both the VLA and the VLBA. 
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5 = +11°17'7.262", with (conservative l-a error of 0.001" in each coordinate). The source is 
unresolved with a size of less than 0.001". 

Ryle Telescope: Observations at 15 GHz with the Ryle Telescope at Cambridge (UK) 
were made by interleaving 15 minute scans of GRB 991216 with short scans of the phase 
calibrator J0509+1011. The flux density scale was tied to observations of 3C48 and 3C286. 

Owens Valley Radio Observatory Interferometer (OVRO): The source was observed for 
a single 13 hr track in two continuum 1 GHz bands (central frequencies 98.481 GHz and 
101.481 GHz) under good 3-mm weather conditions. Gain calibration used the quasar 
0528+134, while observations of Uranus and 3C 454.3 provided the flux density calibration 
scale with an estimated uncertainty of ~ 20%. See Shepherd et al. (1998), for details of the 
calibration and imaging. No source was detected (see Table [l]) . 

James Clark Maxwell Telescope (JCMTQ): Observations in the 350 GHz band were made 
using the Sub- millimeter Common-User Bolometer Array ( [Holland et al. 1999| ). The data 



were taken under good sky conditions on both nights. For flux calibration we used the 
source CRL618, and assumed its flux density to be 4.57 ± 0.21 Jy. The pointing was 
monitored and found to vary by less than 2". See Kulkarni et al. (1999) for details of data 
reduction. The source was not detected at either epoch (see Table [l]). At the position of 
GRB 991216 we derive an average flux of —0.28 ± 1.1 mJy 



3. Results 

In Figure we display the 8.46 GHz light curve, as well as the X-ray and optical 
(R-band) light curves obtained from measurements reported in the GRB Coordinates 
Network (GCN)Q and Halpern et al. (2000). A noise- weighted least squares fit of the form 
F u oc t av was made to each of these light curves. Using all the 8.46 GHz data, including the 
upper limits, we derive a T = —0.82 ± 0.02 (x 2 — 26.5/15; here xl is the reduced x 2 )- 

A similar least squares fit of the optical and X-ray data over the first three days (Figure 
i, yields a Q = -1.33 ± 0.01 (xl = H/28) and a x = -1.61 ± 0.06 (xl = 7.7/3). From a 
more extensive dataset, Halpern et al. (2000) fit a a = —1.071^8 over the same time range. 
We will use their value of a Q in the discussion to follow. The relatively large value of Xr f° r 



7 The JCMT is operated by The Joint Astronomy Centre on behalf of the Particle Physics and Astronomy 
Research Council of the UK, the Netherlands Organization for Scientific Research, and the National Research 
Council of Canada. 

8 http : //lhcawww. gsfc.nasa.gov/docs/gamcosray/lcgr/bacodine/gcn_main. html. 
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fit to the X-ray data presumably reflects the uncertainties inherent in converting the counts 
measured by three different instruments (RXTE-ASM, RXTE-PCA and Chandra ACIS) 
into Jansky flux units. Using the RXTE-PCA data alone avoids this cross-calibration issue 
and yields a x = —1.61 ± 0.05. 



4. The Unusual Nature of the Radio Afterglow: The Failure of the Basic 

Afterglow Model 

The radio afterglow from GRB 991216 is unusual on two counts. First, the radio 
afterglow in the centimeter band does not show the usual rise to a peak value f m (at epoch 
t m ) before undergoing a power law decay. The radio flux appears to decline continuously 
starting from the epoch of the first observation. Thus t m < 1.49 d as compared to the 
10-100 d seen in other bursts (e.g. Frail, Waxman fc Kulkarni 2000| , Frail et al. 1999| , Frail 



et al. 19991) . Second, the temporal decay indices («„) in the radio, optical and X-ray bands 



are markedly different from each other. Proceeding from radio to higher frequencies, a v 
steepens by ~0.4 every four decades in frequency. 

In contrast, the optical and X-ray afterglow appears to find a straightforward 
explanation in the standard afterglow model in which a jet geometry is invoked ([Hafpern et 



\al. 2000| ). Below we show that the radio observations cannot be reconciled with a standard 
jet (or sphere) afterglow model. We then explore possible modifications to the standard 
model. 

The simplest afterglow model is one in which the broad-band afterglow emission arises 
from the forward shock of a relativistic blast wave propagating into a constant density 
medium (|Sari, Piran fc Narayan 1998f) . It is assumed that the electrons in the forward 
shock region are accelerated to a power law distribution for 7 e > j m , dN/dj e cx 7~ p ; here 
7 e is the Lorentz factor of the electrons, p is the power law index and 7 m is the minimum 
Lorentz factor. Gyration of these electrons in strong post-shocked magnetic fields gives 
rise to broad-band afterglow emission. Two modifications to this picture are routinely 
considered. (1) An inhomogeneous circumburst medium (specifically, p(r) oc r -2 ; here p is 
the density at distance r from the source). Such a circumburst medium is expected should 
GRBs originate from massive stars ( |Chevalier fc Li 1999] ). (2) A jet-like geometry for the 
blast wave ( |Sari, Piran fc Hafpern 1999| ). This modification is motivated by the propensity 



of jets in astrophysical sources. 

Regardless of these modifications, the broad band spectrum is composed of three 
characteristic frequencies: u a , the synchrotron self- absorption frequency; u m , the frequency 
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at which the emission peaks (and attributed to the electrons with Lorentz factor 7 m ), 
and v c , the cooling frequency. Electrons radiating photons with frequency > v c cool on 
timescales faster than the age of the blast wave. The evolution of these frequencies is 
determined by the dynamics of the blast wave. The usual ordering of these frequencies at 
epochs relevant to the discussion here is (going from low to high frequencies) u a , v m and u c . 

For GRB 991216 the early radio decay implies that v m is already below the centimeter 
radio band at 1.49 days. The steepening of the afterglow emission from optical to X-ray can 
be explained by placing v c between the optical and X-ray bands. The expected steeping 
Aa is 1/4 which is marginally consistent with a a — a x = 0.54ioio- However, even if we 
ignore this, we are simply unable to explain the decay in the radio band, since no additional 
steepening is expected between u m and v c . 

The standard afterglow model can be made to agree with the light curves by postulating 
an energy slope p which gradually steepens with increasing electron energy j e . We use the 
spherical, constant density afterglow model ( Sari, Piran fc Narayan 1998Q to convert, in 



each band, the observed decay index to p and obtain: p = 2.09 ± 0.03 (radio), p = 2.43io]ii 
(optical), and p = 2.81 ± 0.08 (if v c is below the X-ray band) or p = 3.15 (if v c above the 
X-ray band). We are justified in applying the spherical model for early times (t < t j ~ 2-5 
d) since the jet geometry is manifested only for t > tj ([Hafpern et al. 200C| ). 

Curvature is both observed and modeled in the synchrotron spectra of the non- 
relativistic shocks from supernova remnants which are accelerating cosmic rays (e.g. Baring 
et al. 1999). To date, models of ultra- relativistic shocks favor a universal value of p, 
independent of energy ( |Vietri 2000| , |Gallant et al. 2000]) , but non-linear effects have yet to 
be treated. 

Nonetheless, the invocation of curvature in the energy distribution of the electrons 
cannot explain the observed broad-band spectrum (Figure 0) of the afterglow on December 
18 (corresponding to 1.33 days after the burst). A plausible fit to the entire data is 
obtained with v a = 1.3 GHz, v m = 270 GHz and v c = 7 x 10 16 Hz and f m = 3.4 mJy; this 
fit is displayed by the dashed line in Figure |2|. As the blast wave slows down, v m moves 
to lower values while preserving f m and thus we expect the flux in the centimeter band 
to rise, whereas the observed flux falls. If the afterglow has a jet-like geometry then the 
radio afterglow is expected to rise until the epoch tj, and subsequently decay very slowly 
(/„ oc r 1 / 3 ) until v m passes through the centimeter band, after which we expect to see a 
decline similar to that seen in the optical (/„ oc t~ 2 - 2 ) ([H arrison et al. 1999[) . As can be 



seen from Figure j|, the radio observations are grossly inconsistent with these expectations, 
particularly the decay is much faster than t -1 / 3 . 
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To summarize, while the optical and X-ray observations can be accounted for by a jet 
model, the radio observations are inconsistent with the standard model. This forces us to 
consider afterglow models in which the radio emission (at least in bulk) arises from a source 
other than the usual forward shock. 



5. A Forward and Reverse Shock Model 

The most natural explanation for two components would be an early contribution from 
a reverse shock followed by a forward shock element at later times. This is the explanation 
invoked to account for the early (1-2 day) radio emission from the afterglow of GRB 990123 
( pari fc Piran 19991 |Kulkarni et al. 1999| but see Galama et al. 19991) . The two bursts share 



several common features: in both jet has been deduced with tj ~ few days, both 

were quite bright at gamma-ray energies and finally both had a seemingly small value of t m 
(as measured in the centimeter band). However, in the case of GRB 990123, the peak flux 
of the forward shock was f m < 260 /xJy QKulkarni et al. 1999| ) and the radio light curve was 



dominated by the reverse shock. In contrast, the forward shock for GRB 991216 appears to 
be quite strong. This difference then explains the seemingly different radio light curves. 

At late times (i.e. timescales greater than the duration of the burst) the flux from the 
reverse shock is expected to fall as t~ L8 ( |Kobayashi fc Sari 2000|) . In contrast, the forward 



shock emission rises as t 1//2 for t < tj and then slowly decays, oc t -1 / 3 until the u m moves 
into the centimeter band. Since tj is known from optical observations ( Halpern et al. 2000| ) , 
the remaining unknowns are the strength of the reverse and forward shock emission. 

In this picture, the reverse shock dominates the radio emission for the first few days 
and the model fit consists of mainly fitting a power law with f v oc t~ 1,& . We note that at 
day 1.5, the VLA 8.46 GHz flux and the Ryle 15 GHz flux are comparable. This suggests 
that the reverse shock is already optically thin at 8.46 GHz at this epoch - similar to the 
situation for GRB 990123 ([Kulkarni et al. 1999|) . We deduce the parameters of the forward 
shock by fitting the radio to optical spectrum around t j=5 days to the forward shock model 
(the contribution of the reverse shock is expected to be negligible thanks to the steep decay 
and since t is comparable to tj, the spherical fireball model is still applicable); we find 
!/ ro ~ 1.4 x 10 12 Hz and f m — 1 mJy As can be seen from Figure ^ this reverse-forward 
model provides a reasonable fit to the observations. 

There are two predictions of this model. First, we expect v m to cross the centimeter 
band at i& = tj(u m /8A6 GHz) 1 / 2 ~ 64 d. For t > t b , we expect the radio flux to decline as 
steeply as the optical flux does for t > tj. The low flux values as measured at the VLA 
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around this epoch are in agreement with this model. 

A second prediction is that for t < t&, we expect, the spectrum to rise as v 1 ^ for 
v < 8.46 GHz. Unfortunately, the data are too sparse to rigorously test this expectation. 
Nonetheless, we note that at day 17.44, the spectrum between 1.43 and 8.46 GHz can be 
described by a simple power-law with slope (3 r = —0.45, steeper than a v 1 ^ slope by 3.6a. 
We consider this to be the weakest point of the model but do not consider it fatal since the 
quoted uncertainties include only instrumental errors and do not include external effects 
such as interstellar scintillation. 

The strongest confirmation of this model would have been the detection of an optical 
flash, as in the case of GRB 990123 ([Akerlof et al. 1999|) . The strong radio emission from 



the reverse shock allows us to predict (by scaling from the optical and radio observations 
of GRB 990123, |Sari fc Piran 1999[) that the flash would have been 8th magnitude or even 



brighter. Unfortunately, this event occurred during daytime and therefore was not observed 
by existing prompt optical counterpart experiments (LOTIS - H. S. Park; ROSTE - C. 
Akerlof; pers. comm.) 

We end this section by noting a worrisome and puzzling issue: we are unable to provide 
a consistent explanation for the near-IR, optical and X-ray observations with a standard 
fireball afterglow spectrum. As noted in Figure pj, there is a broad maximum around 2 x 10 14 
Hz - suggesting that this is the peak frequency (u m ) of the fireball. Fitting a template 
afterglow spectrum we obtain the following: u m x = 2.1 ± 0.6 x 10 14 Hz, f m % = 150 ± 10 //Jy 
and u c i = 2 x 10 16 Hz. We note that a similar broad peak in the near IR (and attributed to 
!/ m ~3x 10 14 Hz at At = 0.5 d) was observed for GRB 971214 QRamaprakash et ai. 1998| ). 
However, if we evolve this v m back in time (with v m oc £~ 3-2 ) we predict a rising R-band 
light curve, inconsistent with the observations (Figure [I]). Moving v m to lower frequencies 
solves this problem but we are left with no explanation for the "near-IR" bump. 



6. A Two-Component Forward Shock Model 

We now consider a model in which much of the radio emission arises as the forward 
shock of an additional fireball (hereafter the second fireball). The principal attraction of the 
second fireball is that we no longer need to relate the radio decay rate to those at optical 
and radio frequencies. We clarify that the optical and X-ray observations are explained 
by the forward shock of the fireball (the first fireball) discussed in the previous section. 
As noted earlier, there is good evidence suggesting that the first fireball is a jet. Thus 
the second fireball must be a more isotropic fireball and move at a smaller Lorentz factor. 
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Indeed, in some GRB models, the central engine is expected to inject two fireballs: a high 
r jet and a low T spherical wind. 

A reasonable fit to the radio data of this second fireball (FS 2; see Figure ||]) on day 1.33 
is provided by f m 2 — 1.2 mJy, z/ m2 =7 GHz, and v a2 =2 GHz. The location of the cooling 
frequency z/ c2 is unconstrained. As a test, we evolved the afterglow spectrum forward in 
time. The model does an excellent job reproducing the declining flux density from 1.43 
and 8.46 GHz at 17.44 days (an observation which the reverse-forward shock model fails to 
explain), but at day 60.40 it predicts a 1.43 GHz flux of ~ 100 /iJy, where only an upper 
limit of —57 ± 44 /zJy is measured. Again we consider this 3-cx discrepancy as a major, but 
not fatal, weakness of this model. 

The three inferred parameters (z/ m2 , fm2, v a2 ) allow us to obtain the energy of the 
blast wave and the density of the ambient medium ([Wijers fc Galama 1999|) : E52 ~ 10 2 erg 



and n ~ 10 -4 cm -3 ; these values are relatively insensitive to the value of the unknown v c 
(which is however constrained to lie above the optical band). The large E and small n are 
primarily due to the small value of t m . 

If this interpretation is correct then we have uncovered the first example of a GRB 
exploding in a very low density medium - perhaps the halo of a host galaxy. The dynamics 
of explosions is governed by the ratio E/n, and as noted above, this ratio is perhaps 10 5 
larger than that typically derived in other afterglow. For this reason, both fireballs, the 
high T and the low F fireballs, would then be expanding at high Lorentz factors days after 
the burst. The V for the low T fireball would be an impressive 20 one day after the burst, 
and the fireball would have had a size of 100 microarcseconds three weeks after the burst - 
almost within reach of measurable with VLBI techniques (cf. [laylor et al. 199"7| ). The jet 



fireball would be expanding even faster in which case the opening angle of the jet is not 6° 
( |Halpern et al. 2000| ) but only 1°. 

To summarize, the radio afterglow of GRB 991216 is unusual and cannot be explained 
by the standard forward shock model. A conventional reverse-forward shock model or an 
exotic two-component forward shock model can account for the observations, but each has 
one major (but not necessarily fatal) weakness. Finally, we have no explanation for the 
near-IR bump seen on day 1.33. GRB 991216 shows that there may be yet new surprises in 
GRB afterglows. 
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Table 1. Radio Observations of GRB 991216 



Epoch 


At 


Telescope 




S ± a 


(UT) 


(days) 




(GHz) 


(/*Jy) 


1 QQQ F)pr 1 ft 00 
lyyy i^l^. io.uu 


1 

1 -OO 


nyie 


1 .U 


i i 00+9^.0 

1 lUUXiOU 


1 QQQ Dpp 1 ft 1 


1 AQ 
1 .^y 


VT A 


8 Afi 


qfin+R7 
yuuxu / 


1 QQQ FW 1 8 ^9 
lyyy jjtt. io.oz 




VT RA 


8 AO 


f UOinoO 


1 QQQ Ti&n 1 ft Aft 
iyyy ucc. 10.^0 


1 81 
l.o± 


TPMT 


oou 


OOUznlOOU 


1 QQQ Dor 1 Q ^0 

iyyy ucc iy.ou 


z.oo 


V ) V M w 


QQ Q 

yy .y 


yum / uu 


1 QQQ Dpr 1 G ^ 

iyyy ucc. iy.00 


9 fi8 


VT A 
V l-i 


8 zLfi 




1 QQQ Tic*n 10 4^ 

iyyy ucc. iy.^±o 


9 78 


TPMT 


oou 


— ZUUUxlulU 


1 QQQ FW 90 OQ 
iyyy jjic. zu.uy 


3 49 


rvy le 


10. u 


luuzn^tuu 


1 QQQ Dpr 99 01 
iyyy jjil. ZrZi.ui 


O.O^l 


nyie 


10. u 


— i uznzuu 


1 QQQ T)f±n 9^ ^0 
iyyy l<ci. ZiO.ou 


u.uo 


VT, A 


8 4fi 




1 QQQ npr 94 9Q 
iyyy z^.^y 


7 fi9 


VT, A 


8 4fi 


( mOO 


1999 Dec. 26.40 


9.73 


VLA 


8.46 


170±72 


1999 Dec. 28.24 


11.57 


VLA 


8.46 


211±25 


1999 Dec. 29.43 


12.76 


VLA 


8.46 


136±37 


1999 Dec. 31.26 


14.59 


VLA 


8.46 


123±39 


2000 Jan. 2.01 


16.34 


VLA 


8.46 


130±22 


2000 Jan. 3.11 


17.44 


VLA 


8.46 


131±36 


2000 Jan. 3.11 


17.44 


VLA 


4.86 


126±31 


2000 Jan. 3.11 


17.44 


VLA 


1.43 


257±51 


2000 Jan. 6.15 


20.48 


VLA 


8.46 


123±30 


2000 Jan. 23.95 


38.28 


VLA 


8.46 


79±31 


2000 Jan. 28.16 


42.49 


VLA 


8.46 


148±33 


2000 Feb. 05.18 


50.51 


VLA 


8.46 


3.1±30 


2000 Feb. 15.07 


60.40 


VLA 


1.43 


-57±44 


2000 Feb. 15.07 


60.40 


VLA 


8.46 


9.6±24 


2000 Mar. 03.85 


78.18 


VLA 


8.46 


47.0±19 



Note. — The columns are (left to right), (1) UT date of the start of each observation, (2) time elapsed 
since the GRB 991216 event (i.e. i = 1999 December 16.67 UT), (3) telescope name, (4) observing frequency, 
(5) flux density of the radio transient, with the error given as the rms noise on the image. The epoch on 
Jan. 23.95 UT is an average of two days of data (Jan. 21.95 UT and Jan. 25.94 UT). All VLA observations 
were made when it was in its "B-array" configuration. 
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Fig. 1. — Broad-band light curves of GRB 991216. Upper limits are plotted as the peak 
flux density at the location of the afterglow plus two times the rms noise in the image. The 
R-band data are taken from Halpern et al. (2000). Optical magnitudes were converted to 
Jansky flux units ([Fukugita, Shimasaku fc Ichikawa 1995|) but no correction has been made 
for Galactic extinction. The X-ray data are measurements taken by the ASM (*) and PCA 
(filled squares) instruments on RXTE ( Corbet fc Smith 1999| , PTakeshima et al. 1999|) , and 
the Chandra X-ray Observatory (□; Piro et al. 1999| ). X-ray fluxes are converted to Janskys 
using the X-ray slope j3 x = —1.1 derived by Takeshima et al. (1999). The solid lines are 
noise-weighted least squares fits to the data, with the slopes a u as indicated (see text for 
details). 
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Fig. 2.— Radio to X-ray spectral flux distribution of GRB 991216 on 1999 December 18.00 
(Ai=1.33 days after the burst). Optical and infrared measurements are taken from Halpern 
et al. (2000) and are interpolated to this epoch assuming the decay rate of a Q = —1.07^11, 
as measured by Halpern et al. at early times. The optical/IR data have been corrected for 
Galactic foreground extinction flSchlegel, Finkbeiner fc Davis 1998|) , giving E(B— V)=0.634 
with an an uncertainty of 10%. The 1.4 GHz upper limit (plotted as three times the rms 
noise) and the 4.8 GHz data point were taken at the Westerbork Synthesis Radio Telescope 
(WSRT) by Rol et al. (1999). The flux density at 8.46 GHz derived by extrapolating of 
the power-law decay in Figure [I]. The upper limits at 100 GHz and 350 GHz have been 
extrapolated back to this epoch by assuming a worst case decay rate of a Q — 2a. The dotted 
and solid lines are fits to the data for a synchrotron spectrum from a relativistic blast wave 
as specified by Granot et al. (1999a; see their Figure 10 for the equipartion field model). We 
assume p = 2.2 and scale it by v m and f m to derive a function g(v) with asymptotic limits 
of z/ 1 / 3 and z/ p_1 )/ 2 . We account for synchrotron self-absorption at v a by multiplying g{y) 
by F„ = [1 — exp(— t)]/t, where r = (z//V a )~ 5 / 3 ( Granot, Piran fc Sari 1999b]) . 
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Fig. 3. — Observed and model light curves at 8.46 GHz. The dot-dashed line is the power- 
law fit from Figure |]. The thick solid line is the two component model discussed in §|| 
consisting of a reverse shock (dashed line) and a forward shock (thin solid line). See text for 
more details. 



